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' The effects of including compressibility in the Kemp-Sears 

problem of aerodynamic interference between moving blade rows are 
examined. Methods of linearized, subsonic, plane, unsteady flow are 
adopted. The major new effect is that a resonance appears at certain 
combinations of flow Mach number, tip Mach number, and blade 
vane ratios. The resonance is at exactly the Tyler-Sofrin cutolT condi- 
tion for rotor-stator interaction. At such conditions the unsteady lift 
on a blade row due to externally imposed nonstationary upwash 
vanishes. However, the resonance appears to be very sharp and 
seems to be more significant as an indication that annind this condi- 
tion the unsteady lift changes very rapidly. 


This paper is concerned with unsteady blade forces developed on an 
axial compressor blade row due to unsteady periodic wave disturbances. 
The applications discussed pertain to unsteady disturbances produced 
on a blad(^ row due to steady-state lifts of adjacent rows moving past the 
row of interest and due to viscous wakes sht'd by upstream rows. 

The framework for the analytical approach to the solution of this com- 
plex problem is contained in two pioneering papers by Kemp and Sears 
(refs. 1, 2). Kemp and Sears adopt the n-presentation of blade wheels as 
infinite cascades of two-dimensional airfoils. They consider a typical 
airfoil in the blade row on which unsteady forces are to be calculated as 
an isolated airfoil in linearized, unsteady, incompressible flow. This 
isolated airfoil is subject to unsteady disturbances from various sources. 
Reference 1 considers nonstationary upwash due to translation of steady- 
state design fields of an adjacent row. In this case the effect of a row upon 
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an adjac(‘nt row downstream of it and the nonstationary upwash due to 
the “inviscid” wakes is included, lleference 2 considers nonstationary 
upwash due to viscous wakc's shed by an upstream row. In tliis paper two 
refinements are considered, k'irst, a typical airfoil of the row on which we 
wish to calculate the unsteady blade forces is regarded as a member of an 
isolated, infinite cascade of airfoils rather than as an isolated airfoil. 
Second, in calculating botli the nonstationary upwa.sh and the un.steady 
forces we include compressibility effects. INIethods for linearized, subsonic, 
plane', compre'ssible un.stc'ady flow are' adeejite'd. For the jire'sent, we have 
not included neinstationary upwashes due tee “inviscid” wakes. The 
motivatiein feir ince)r])e)rating the'see refinements is simply that in meest 
present-day applications whe're' unsteady blade force's in fans and com- 
jeresseers ;ire‘ eif inte're'st, the Mae'h numbe'rs e>f the' flow are too high for 
ceempressibility effects to be' negligible'. 

('onside'r an isolate'd, infinite flat-plate cascade' of ide'iitical airfoils 
space'd s ajiart wheese cheerd lengths havee tie'e'ii normalize'd tee unity and 
with flow at Mach numbe'r M through the cascaele. M is restricted te> 
M <\ (subsonic flow through the- blade' ])a.ssage'.s). d'lu' {irobh'm i.s to 
calculate the un.ste'ady lift on a typical blade' e>f such a e'ascade due to a 
known nonstatieuiary upwa.sh een it whe're' the' tiine'-eh'penele'nce' e>f the 
neeiYStationary upwash is e>f simiih' harmeenic ty]>e'. The' nem.statieinary 
upwash will later be re'late'el tei the' transliitie)ii e>f aeljace'ut ste'ady-state' 
de'sign fie'lels anel to vise'eeus wake's she'd by an ujistre'emi reew. Due to the 
symmetry eif the cascaele' the' neenstatieenary u])wash een the' ath blade is 
esse'utially the' same' as ein the' zeroth blaele, e'xe'ept feer a ]fliase differe'uce 
facteir of exp {jay)', wlu'ree we; will late'r re'late' y to the aerodynamic 
parame'ters. 

The me'thoel e>f seelution adeeptc'el is the me'theed of distributed singulari- 
tie's and singular inte'gral e'eiuatieens. An x — y ceeeerdinate system is use'd 
as sheewn in figure' 1. Le't ^ de'nefle' the' running coeerdinate on the' zereeth 
blade that runs from —5 to Le't the' unste'ady lift di.stribution on the 
zeroth cheerd be de'iieite'd by Le't A',.(x,^) deneite a kerne'l functieen 

that gives the nemsteitieenary upwash at the' peiint (a",0) een the' ze'roth 
blade due tee an infinite' row e)f e'e[ual oscillating feerces e>f unit .strength with 
phase shift y located at | and its corre'sponding points; i.e., at 

^+ 7 is sin (os), as cos {a,) 

whe're m = 0 , ±1, ±2, . . . , e'te'. is known since it is merely the 

sum (>{ upwash contributiems due tej unit force's of kneewn phase and 
locatiein. Led Vd(x)e^‘^‘ be the kneiwn nonstatiemary upwash due to the 
adjacent rows. Then the boundary ceendition that the ve'locity normal to 
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the blade chords hi' zero leads to the integral e(|iiation: 


f 

J _i /o 




(1) 


where K^{x,^) and i\i{x) are known. Only F{^) is unknown, and, after 
solving equation (I), the total unsteady lift is found from 


(total lift) 


An 

= / 

*'- 1/2 


(?) 


( 2 ) 


A',.(x,^) has a singularity; i.e., a term as \ /{x — ^) \ hence, equation (I) 
constituti's a singular integral eiiuation for the forc(' distribution. There 
are two other reciuinanents on F(^) \ first, that it vanish at the trailing 
edge ^=5 (the Kutta-Joukowski condition), and, second, that it have a 
square-root singularity at the leading edge; i.e., the force distribution 
tends to =» as 5 as l/\/ 


CALCULATION OF KERNEL FUNCTION 

The physical significance of the kernel function K,,{x,^) is repeated: 
It is the upwa.sh at point (x,0) on the zeroth blade due to an infinite row 
of unit forces located at 

^-\-ns sin (as), ns cos (as) 

with n = 0, ±1, ±2, . . . , etc., where the jihase of the oscillating force 
on the /ith blade is related to that on the zeroth blade by a factor 
exp ( jny) . 


516 UNSTEADY FLOW AND NOISE 


In the derivation of the kernel function we have relied heavily on a 
recent paper by Kaji and Okazaki, who treat an analogous problem 
(ref. 3) . We are espi'cially indebti'd to their repeated and ingenious use of 
the Poisson summation formula to switch from infinite series of Hankel 
functions to infinite series of exponential functions. In view of their paper, 
we omit manj^ of the details of the derivation of the kernel function and 
follow Kaji’s notation. 

The first stc'p is to note that the nonstationary upwash Vn{x,0,t) at 
(x,0) due to a force of comph'x strength exp ( jny) located at 
Kf+?is sin (as), ns cos (as)] is as follows. 

Let k = u/a, ^ = \/ l~M-, p = density of uniformly flowing medium, 
a = speed of sound, a;,. = x — ^— ns sin (as), ?/„= — ns cos (a»), and 

Xn' = x' — ns cos (a,). Then 


F„(a;,0,0 I 


kH 


V l-M- 


. kMx„ 


X 




^ Vxn^+^-y,A 
^ j 
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L/3-a 


//o«) 


k^ exp { ~j — Xn 




and 


/ ^n 

exp 


■ . 

71/(l-il/-’) 


— y/xn'^ + ^hjn- dx 
fi-a 


(3)'^ 


Kv{x,^)= X Vn{x,0,t) 


( 4 ) 


The above expression for V„ may be dt'duced from ecpiations (8) and (11), 
section 14.3, of Y. C. Fung’s “An Introduction to the Theory of Aero- 
elasticity.” Thus, 


and are Ilankel functions of the second kind of orders 0 and 1, 

respectively. 
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Kr(X ,0 = 
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(5a) 


7 = y (2i m+P)/ //2) 
'“i V(2«i+P)^+r 


[tc-s/ (27m+P)2+r;2] exp [j(2im+P)5] 


(5b) 


h= S Po®[wV (27m+P)2+)j2] exp [](27m+P)5] 


— X 


(5c) 
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li represents the contribution from the integral term in equation (3) ; 
we give below the transformed form of it after application of the Poisson 
summation formula; 


h = iY, 


exp \_jnP — t )\/ (5 — «)- — 


-\/ (5 — /i)-— w- 


j{n-b)P+j 


k 


M(l-JP) 




for (x — ^) < 0 . 


(5d) 




{5 — n)- — iv- 


exp [_jnP — r ]\/ {8 — n)- — w-^ 
j{n-5)P+j 1/(^1 IP) 


2 r] exp 




k{x-0 


(n — S)Pd 


k 

.1/(1 -dP) 


+ jj-[( (5 — — iP] 


for {x~i) >0. 

For X 9 ^^, the expressions for I\, 1 2 , and 1 3 may be transformed l)v 
Poisson’s summation formula as under 

7i = -^ X exp — — tP] (5«) 

Tll'T] 

h = — £ CjnP-vV(S-fi)'^-P-J (5b) 

TTlC V (5— 

J «= I 

/3 = - 2] - 7-= - exp [jhP-jj (5c) 

The expressions for P, P, and I 3 diverge as x~^^ in equations (5a), (55), 
and (5c) , or as P, ij— ^0, which is merely the indication of the 1/| x — ^ \ and 
log ] x — ^ I types of singulariti(\s in K,.{x,^). To locate these singularities, 
we examine the small value expansions of the Hankel functions as follows. 

1-— (70- log 2)1 


//o®(^) ~ — “ log (a;) + 

7T 


where 70 = 0.5772 = Euler’s constant and (x) ^2j/Trx. With the aid 
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of the above it is easily shown that 


K,{x,^) ■ 


jk 


2j/3* 


4TTp(3-Ma \_k{x~ M 
is regular as x-^^. This motivates writing 


log (I x~^ I) 


= Kr{x,^) ( 6 ) 


/3 k log (I x — ^ 

K,.{x,0 = K Ax, - i-Tp- 

2TTpMa (x* — ZTTpfiM-^a 


Finally, let 


-2pMa 

fix) = i’d{x) 
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Equation (1) noiv reads as 

F(^) 


•'-m ^ 


jTlog (1 i) 1 


2pM (ITT 

^ ' M ( f- AT) ^ “ Jx^ 


(IA-=f(x) (9) 


CALCULATION OF UPWASH DUE TO TRANSLATION OF 
ADJACENT STEADY-STATE FIELDS 


Potential Interaction 

The probk'm of potential interaction (see fig. 2) is to calculate the 
nonstationary upwash on line AB due to steady-state design lifts on a 
blade row on the right-hand side, which translates downward at a ei'rtain 
speed corresponding to a tip Mach numlxT M,. In a linearizi'd tri'atment, 
clearly the stagger {a^ and Alach number Mr of flow through the ad- 
jacent row are related to M,, M, and as by 



row. 
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«r= tan ' 


Mt + M sin {a,) 
M cos («s) 


and 

Mr= Vd/r+dr- + 2d/,d/ sin (a.) 


Let Trr = rtd/r. W(‘ first find the solution for velocity comjioiH'iits u' and 
v' , parallel to the x'-y' coonlinate system of figure 2, du(‘ to ecpially sjiaced 
concentrated unit force's at the origin and its corresponding jioints as 
shown in figure 2. (Note' that tlu' blade exe-rts a force' een the' fluid e'epial 
and eipposite' fee the' lorce' e'xe'rte'd by the' fluid een the' blaele'.) We' use a 
frame' eif re'fe're'uce' fixe'el with re'spe'ct te) tlie> translating blaele' reew se> that 
we' have' a ste'ady-state' preebh'in. We' heive' te> e'eenside'r the' effe'ct eef a 
sum e)f force's; 


1 X) 

X \ 


where 1 denotes a unit feerce' vector anel 5 stanels fe>r the' Dirac ele'lta 
function. By using a re'sult eai page (i8 eef re'fe're'tice' 4 ceince'rning the sum 
of an infinite renv of eepially sjeace'd de'lta functions, we' finel that the 
above is clearly equal to 



1 + 2 X''XP 

7i = l 


270 )?/^ 


since our use e>f comph'x Ibrms always im])lie's that re'al jiiirts are' te> be 
taken. The linearize'd eeiuatiems eef meetieen and ceentinuity are' 
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(10c) 
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In what follows wc omit the %j' indepemdont t(>rm in the force since it gives 
a stationary upwash (cf. p. 591 of ref. 1). In equations (12a), (126), and 
(12c), p' and p' stand for small perturbations of the demsity and pressure. 
Eliminating // from equations (106) and (10c) we derive 


cos (or) 


-- + sin (a,) 
dx 
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00 
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= 0 


Since u' , v', and the d(dta function term vanish far from the blade row. 
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(E utta-Joukowski Law). 

Next we eliminate the force terms in equations (126) and (12c); 
assuming an isentropic relation between p' and p', we obtain: 

[1-dE^ cos= («r)]+~ [l-3E^si.P (a.)] 

dx dy 


= d/r-sin («r) cos {Ur) 


'dll' 

Jy' 


d_i^' 

dx' 


(12) 


(modifu'd continuity equation). 

Using equations (11) and (12), single equations for u' and v' may be 
obtained that may be solved by recpiiring that u' and v' vanish as oo . 

We omit the details and give the result. 


u' for x'^0 = 
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( 136 ) 


The upwash normal to line AB (fig. 2) is 
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Vd = v' cos (as) —u' sin (as) (13c) 

Thus I'd for x' ^ 0 is 

=F1 r sin (ar) cos jar) sin (as) j sin (as) — 

prffi'r 1 — jT/r- COS^ (ar) ^ 1 — COS^ (or) 


00 

xX exp 
1 



2 irnx' sin (a,) cos (ar) =F -\/l — jl/r^ 
d 1— 3/r^ cos- (ar) 


(13d) 

The effect of distributed loading on a finite chord may be estimated by 
integrating the results of equation (13d) over the finite chord. We again 
omit the details since the derivation is very similar to that outlined on 
pages 589 through 590 of reference 1. Noting that these upwash fields 
translate with respect to the adjacent rows, om? readily obtains the non- 
stationary upwash on the adjacent rows. 


Viscous \\ ake Interaction 

The form of the nonstationary upwa.sh contributi'd by the viscous 
wake intraction mechanism was assumed to be the .sam(> as in ri'ferenci'- 
2. The pertinent eiiuation giving the upwash is ecpiation (28) of nderenco 
2 . 


FINAL FOR.MLLATION OF INTEGRAL EQUATION (9) AND 
METHOD OF SOLUTION 

From figure 3, clearly the frequency co of unsteady lift is {2vn/d)M, a 
and thus k = u/a= {2wn/d)M In what follows we consider (>ach harmonic 
n separately. The phase lag y is discussed on pag(> 592 of reference 1. It 
is easily shown that 

Assume for F{^) the form 

© N 

+ 2] sin (n0) 

1 

where ^ = — 5 cos (^, 0 = 0 at the leading edge, and </> = tt at the trailing edge. 

All the above terms are zero at the trailing edge (Kutta condition) , and 
the series has the usual square-root singularity at the leading edge. 
Ao, Ai, A 2 . . . A„ are, of course, unknown. We let x = — 5 cos 6 so that 
6 = 0 corresponds to x = — | and d = ir corresponds to x = | and denote by 
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UPSTREAM 

ROW 


DOWNSTREAM 

ROW 



MfO 


SOLIDITY CTu, 
STEADY-STATE 
LIFT Fu 


SOLIDITY CTji, 
STEADY-STATE 
LIFT Fd 


Figukk 3 . — Sketch for potential interaction due to two adjacent rows. 


G(n,d)= I COS {n<t>)Kr{x,^) (Itj) 


whore ?i. = 0, 1, . . . , A' + l. 

The G(n,d) will he evaluated numerically by Simpson’s rule and will 
pn'sc'iit no difficultic'S since Kr{x,i,) is houndc'd. Then it is easily shown by 
using th(> nvsults of equations (16), (17), and (18) of the appendix that 
equation (9) may be written as 
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The above equation holds for O<0<7 t, and the method of collocation of 
points involves satisfying the above equation exactly for (A’ + l) equally 
spaced values of 6 between 0 and tt and thus deriving (A’ + l) simul- 
taneous equations for Ao, Ai, , A.v, which may be solved by matrix 
inversion. 

The net lift is 



and the magnitude of unsteady lift is the amplitude of Tr/2[Ao-|- (Ai/2)]. 


NUAIERICAL DETAILS AND TYPICAL RESULTS 

The unsteady lift distribution was assumed to havc^ the form of an 
term series. I'lic check on whether the chosen 4/ is satisfactory 
is whether the terms Ao, Ai, . . . , resulting from the solution of the 
simultaneous equations (14) converge rapidly enough. Judging from the 
calculations performed in this papc'r, the value of M to be used increases 
with the Alach number of the flow through the blade passages. Up to a 
Alach numb(>r of about 0.5, M = 7 suffices. Between Mach numbers of 0.5 
and 0.8, J/ = 11 sufficivs. Beyond a Mach number of 0.8 it seems necc'ssary 
to use M= 15 to get good convergence. The use of the present analy.sis for 
Mach numbers close to unity is not very valid anyway because for such 
high IMach numbers the convected wave equation (from which eq. (3) is 
derived) is not a valid linearized equation for describing the nonsteady 
flow. 

In figure 4 we have plotted results for potential interaction'* on a row 
with flow at a IMach number of 0.1 due to a row downstream. This case 
should be analogous (owing to the low IMach numbers) to a case cal- 
culated in figure 5 of reference 1. The Kemp-Sears results and results of 
this paper compare reasonably well. 

A similar check with the Kemp-Sears results is obviou.sly de.sirable for 
viscous wake interaction. In reference 2, in the interest of obtaining a 
closed-form solution, the upwash used to calculate the unsteady lift is 
taken at selected points on the airfoil. Two sets of results pertaining to a 
stator rotor sequence as sketched in figure 5, one corresponding to the 
upwash at quarter chord from the k'ading edge and another corn'sponding 
to the upwash at quarter chord from the trailing edge, are presented in 
table 1 of reference 2. The methods used in this paper make such an 


'■For all potential interaction calculations reported in this paper (Figs. 4, Go, 7h) 
the steady lift distribution is assumed to be of the flat-plate tj'pe (see p. .')94 of ref. 1). 
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Figure 4. — Potential interaction due to downstream row. Comparison with Kemp-Sears 

(reference 1). 



Figurk o. — Stator rotor sequence ttsed to compare 
results for viscous wake interaction with Kemp- 
Sears (reference 2). 


approximation uiiiipcossary ; our results are compared with the Kemp- 
Sears results in table I below. The calculations in this paper are again for a 
row with flow at a Mach number of 0.1. 

From table I it is seen that the first harmonic results compare well if 
Kemp-Sears results corresponding to an upwash at a quarter chord from 
the leading edge are used. 
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In figure 6a we have plotted results for unstcadj' lift assuming the 
velocity triangles of figure 4 to be linearly scaled up in i\Iach number. 
This could be repre.sentative of changes in unsteady lift as one runs up a 
fan or compressor on a constant loadline. 

The forms for 1^, I 3 , and /j obtaiiu'd b 3 ’ the use of the Poisson summa- 
tion formula (cqs. 56, 5c, and 5d) indicate that if, for anj' integer m, 

(6-m)2-o)2 = 0 . . . (14) 

then h, h, and /j— ^= 0 . Since in equation (1) involves h, h, and 

li, this means that if ecjuation ( 14) is .satisfied then X,.(x,?)— > 00 . The only 
way in which one can obtain a bounded v,i{x) under the condition that 
A'v(x,^)— >03 in (>quation (1) is to have F{^)—*0. Thus the resonance 
condition denoted by ('quation (14) is one for which the unsteady blade 
forces vanish. Physically this condition ari.ses when purely transverse 
waves are produced in the blade jia.ssages (i.e., waves traveling only in the 
tangential direction). At this resonance condition, waves emitted from 
one surface, say the upper surface of a blade, travel transversely and 
arrive at the lower surface of an adjacent blade with the time of travel 
being such that the incident wave phase is exactlj' antiphase (i.e., with a 
phase difference of some odd intc'ger multiple of 180°) with the phase of 


Tahle I. — Viscous Wake Interaction' 


Keinp-Sears results 


Results of jiresent pajicr 


Harmonic 

number 

Co 

X^p\va.sh 
evaluated at 

Co 

Harmonic 

number 

Cd 







1 

0.01 

Q.C. from 
L.E 

0.029 

1 

0.01 

Q.C. from 

T.E 

0.023 


0.02 

Q.C. from 

L.E 

0.0.')" 


0.02 

Q.C. from 

T.E 

0.04.5 


Q.C. = quarter chord, ('n = profile drag coefficient, L.E. = leadiiis edse, T.E. = 
trailing edge, C/. = coefficient of unsteady lift. 


> Comparison with Kemp-Sears results for case sketched in figure 5. 
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Figure (ia. — Potential interaction for di/namicalli/ similar velocity triangles. Effect of 

Mach number. 


the sources on the lower surface of this adjacent blade. This phase can- 
cellation prevents development of any unsteady lift. 

The resonance condition of ecjuation (14) may be rewritten in terms of 
tip !Mach numbers and flow Mach numbers as 


M,= 



co.s^ (a») — fl/ sin (cis)] 


(15) 


For velocity triangles of the type shown in figure 4, ilf < = M \/2, a* = — 45°, 
and the solution of (15) yii'lds resonant Alach numbers of 


M = 


(?w-l ) 
\/ wF-fl 


where w is a positive integer. The first significant resonant ^ilach number 
M is thus 


M = —7= = 0.447 
V5 


and the second is 
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These two ordinate lines (corn'spoiidinj;; to il/ = 0.447 and 0.G32) are 
shown ill figure' (ie as first and se'cond re'sonanee' jioints. As can bo noted, 
th(' resonanci' is oxtri'mely shar]) in that tlie ajijiroach of the unsti'ady lift 
to zero as 4/— ^1/^/5 is extremi'ly sharp, being ri'pn'si'nted by an almost 
vi'i'tical droj) in figure G«. 

The results of figure (in suggi'st that the re'sonanci', while undoubtedly 
denoting a point of zero unsti'ady lift, is much too sharp to have practical 
significance as a condition of low unsteady lift. Jlowever, the resonance 
points do have' ceinside'rabh' significaiu'e' (as may be- eibse'rve'd from fig. (in) 
as de'liiH'ating rathe'r diffe're'iit familie's of variations of unsie'ady lift with 
■Mach number. 

A similar result is shown by Kaji and Okazaki in refere'iice 3. They 
ceinside'r in refeTe'ne'i' 3 (he effe'e'f of a flat-j)late' cascade' with flow on an 
incide'iit sounel wave' as sheiwn in figure' (ih. An unste'ady force' distributiem 
on the blades tei cance'l the vi'leie'itie's ineluce'd by the' incide'iit sound wave' 
is sought as the solution tei an inte'gral ('epiation eif the same type as 
equatiein (1). Late'r the' e'ffect of the' unste'ady force' di.stributions is 
integrated tei eibtain far-fie'ld pressure (sound) waves in the transmitted 
and reflected re'gions. The- cascade in figure' Gh behaves as a diffraction 
grating with respect tei the incieh'iit senind wave. One iireipagating rc'flected 
and one preijiagating transmitted wave arc always produce'd by the- inte'r- 
action of the incide'iit sound wave and the'iilade row. The' reflected wave 
corresponds to a specular reflection of the incident wave by the blade row. 



Fku'HF, (il). — Configuration of cascade iciiosc results arc shown in figure 6c, of 

reference 3. 
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and the transmitted wave has the same orientation as the incident wave. 
The basic transmitted wave, produced under all circumstances, is denoted 
by the authors of riderence 3 as a r = 0 mode^ in figure 5 of their paper (the 
upper half of which is reproduced as fig. 6c in this paper). Higher order 
modes arc also produced if, for the interaction of the incident wave and 
the blade row, more than one resulting mode is above “cutoff." In figure 
6c, for a less than about 40°, a basic transmitted mode (labeled r = 0) 
and a higher order mode (labeled i/=— 1) are produced. For a greater 
than 40°, the basic transmitted mode and a higher order mode (labeled 
r=+l) are produced. The pressure transmi.s.sion coefficient is the ratio 
of amplitude of transmitted wave to the amplitude of incident wave. The 
ordinate around a = 40° in figure 6c represemts a resonant condition of the 
type of equation (14). In the example of figure 6c, the orientation of the 
higher ordiT transmitted wave undergoes an abrupt change* as one passes 
across the resonant incidence angle. Note how the pressure transmission 


PRESSURE TRANSMISSION COEFFICIENT FOR 
CURVE I:V =0 MODE 
CURVE n- F = I MODE 
CURVE ED- F =-l MODE 



® Corresponds to m — 1 = 0 in this paper. 
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coefficient associated with the basic transmitted wave ((/ = 0 mode) under- 
goes an abrupt change as one crosses over the resonant incidence angle. 

This type of resonance' effect was apparently first observed by Runyan, 
Woolston, and Rainey (ref. 5). Their concern was with the effect of wind- 
tunnel walls on the lift forces developed by an oscillating wing in two- 
dimensional subsonic compressible flow. The resonance phenomenon was 
experimentally confirmed by these authors. 

In view of the very rapid change of events just around resonance, it 
might appear worthwhile to attempt an analytical solution of equation 
(1) close to resonance. However, McCune, in a different context, has 
examined the merits of replacing cylindrical wave functions (which should 
be employed in a proper three-dimensional analysis) by two-dimensional 
approximation in such problem.s (ref. G). He .shows dearly tliat such an 
approximation breaks down precisely at these resonance frequencies. 
Thus, it seems of dubious advantage to pur.sue the cascade jdane analysis 
any further near rc'sonance. Finally, one may easily show that the reso- 
nance conditions are precisely the conditions at which succe.ssive rotor- 
stator interaction modes of the classical Tyler-Sofrin analysis (ref. G) 
are cut off. Thus the cutoff frequencies introduced by Tyler and Sofrin 
as delineating regimes of acoustic propagation or decaj’ of successive 
interaction modes are sec'n to play an important role in the estimation of 
the unsteady blade force problem. 

From equation (13d), one notices that the; exponential decay rate of 
the potential flow field of an adjacent row is altered from its incomprc'ssible 
value by the factor 

yr^M? 

1—M/ COS^ (ttr) 

To show the effect of this factor, we have' jilottcd the decay of unsteady 
lift with spacing for o:r = 0, 45°, and 60° in hgure 7. 



UNSTEADY LIFT, 1st HARMONIC 
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ALL DOWNSTREAM ROWS, STEADY-STATE LIFT=I.OO 



( I ) DOWNSTREAM ROW 
STAGGER = 0° 


( i i ) DOWNSTREAM ROW 
STAGGER = 45° 


(iii) DOWNSTREAM 
ROW STAGGER 


Ficcuk 7a. — Effects of stagger angle on decay of potential flow fields — configurations 

considered. 



Fioure 7b. — Effects of stagger angle on de- 
cay of potential flow fields — results for 
unsteady lift for configurations of figure 
7a. 


0.01 
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APPENDIX 


Two well-known principal value integrals are 


and 


1 r 

COS (?((/)) dip 

T Jo 

_cos (0) — cos (6) _ 


1 /■’" sin {ii4>) sin f</>) (l4> 
IT 2|) cos (9) — c((s (4>) 


sin {n9) 
sin 9 


cos {u9) . . . 


( 16 ) 


(17) 


An expansion ofjk (log | x — ^ |)/2.1/( I — .1/-) convergc'iit ior 97^4>\a 


-J k 

M(l-M-) 


^ cos {ti9) cos 

iog(2)+i: - — 


(18) 
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DISCUSSION 


A. AHDELHAMID (Carleton University); I would like to eommiTit 
on the effi'ct of the blade being in a row or being isolati'd and on whether 
it’s a first- or a second-order effect. I believe it’s a si'cond-order effect 
because the steady vorticity on the blades of the same row as the blade 
that you are considering will not contribute to the contour it imposes. 

J. E. FFOWCS WILLIAMS (Imperial C'ollege of Science): It seems to 
me that when you have this guaranteed repetitive .system you must 
gimerate a re.sonance. You’ve got to pick them up .somehow. My point of 
concern is the way in which the boundary conditions are put in. There are 
difficulties concerned in judgments, and these turn up when the re.sonant 
wave fronts are going at a parallel, one normally observes. Now if thiw’re 
going parallel to the surface, the procedure to adopt is one to control the 
velocity on the surface. If one asks what pressure is reipiired on that 
surface to bring about a control on the velocity, it turns out to be infinite 
by the way his blade is set, so a more realistic boundary condition for 
any practical system would be a jiressure-releasc' condition. 

MANI (author) : Two points have been raised in the discu.ssion period. 
The first concerns the proper boundary condition to be usial near the 
resonant (cutoff) frequency. The author agrc(>s that if very large pres- 
sun's result the assumption of perfectly rigid blades is not a suitable one 
and should be replaced pi-rhaps by an impedance condition. It is worth 
reiterating, however, that from the point of view of the isolated two- 
dimensional cascade model’s n'presentativeness of the actual situation 
in the turbomachine, it would not be worthwhile to pursue matters much 
further in the cascade plane near cutoff. The effect of adjacent blade 
rows, open-end terminations, and three-dimensional effects become all- 
important at this condition. 

The second pertains to the effect of adjacent blade rows. It should 
be noted that estimates of effect of adjacent rows as carried out in the 
original Kemp-Sears papers is not entirely sufficient. Kemp and Sears 
assumed the fluid to be incompressible and, with this assumption, all 
iionaxisymmetric flow patterns exhibit exponential axial decay. When 
the fluid is regarded as compressible, some of these flow patterns (those 
above cutoff) exhibit no axial decay and hence the treatment of the 
problem on the ba.sis of isolated blade rows subject to given external 
upwash becomes questionable. 


